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Abstract 

We study the possibility of identifying dark matter properties from XENON-like 100 
kg experiments and the GLAST satellite mission. We show that whereas direct detection 
experiments will probe efficiently light WIMPs, given a positive detection (at the 10% level 
for m x < 50 GeV), GLAST will be able to confirm and even increase the precision in the 
case of a NFW profile, for a WIMP-nucleon cross-section cr-^—p ^ 10 ^ pb. We also predict 
the production rate of a WIMP in the next generation of colliders (ILC), and compare their 
sensitivity to the WIMP mass with the XENON and GLAST projects. 



1 



Contents 

1 Introduction 

2 Direct detection 4 

2.1 Differential event rate 4 

2.2 The statistical method 5 

2.3 The XENON experiment 6 

2.4 Influence of astrophysical/background assumptions |8 

3 Indirect detection \w 

3.1 Differential event rate 10 

3.2 Modeling the galactic center background 11 

3.3 The GLAST experiment 12 

3.4 Varying background parameters 13 

3.5 Impact of Final States \16 

4 Direct versus indirect detection experiments LL7 

5 WIMPs at present and future colliders |^ 

5.1 The Approach 20 

5.2 Basic Results 22 

5.2.1 Non-polarized beams 22 

5.2.2 Polarized beams [25 

6 Complementarity [25 

7 Conclusions 26 



2 



1 Introduction 



There exists strong evidence that a large fraction of the matter in our Universe is non-luminous 
PQ. Such evidence includes the motion of cluster member galaxies [2j, gravitational lensing |3j, 
cosmic microwave background [4], observations of the flat rotation curves of galaxies [5], etc. 
Dark matter plays a central role in current structure formation theories, and its microscopic 
properties have significant impact on the spatial distribution of mass, galaxies and clusters. 
Unraveling the nature of dark matter is therefore of critical importance. A Weakly Interacting 
Massive Particle (WIMP), with mass lying from the GeV to the TeV scale, is one of the 
preferred candidates for the dark matter of the Universe. 

Different experimental programs are developing huge efforts to observe and identify the 
particle nature of dark matter. This can be achieved by direct measurement of the recoil 
energy of a nucleus when scattered by a WIMP, or indirectly via the observation of WIMP 
annihilation products. In both cases, the sensitivity depends strongly on the background and 
on the theoretical assumptions of the model. It would be interesting to combine all these efforts 
to invent intelligent strategies for determining the nature of dark matter [6J. Recently, several 
works (see, for example, Refs. [7] and [8] for the case of direct detection and [H [10] for the 
indirect detection case) have shown that precision measurements of the mass of WIMPs are not 
only reserved to the domain of accelerator physics. In all of these studies, model independent 
bounds are derived for annihilation cross-sections, masses or WIMP-nucleus scattering cross- 
sections. The drawback of a model-independent framework (lack of determined microscopic 
processes) is largely compensated by the universality of the method: instead of restricting 
a theoretical parameter space, we restrict observable physical quantities (masses, branching 
ratios). Indeed, these limits are valid for a great number of WIMP candidates such as, for 
example, the supersymmetric neutralino, the lightest Kaluza Klein excitation, etc. 

The aim of the present work is to analyse two of the most promising experiments, XENON 
[TT] and GLAST [12], calculating and comparing their sensitivity to a WIMP mass depending 
on the astrophysical hypothesis (velocity distribution of WIMPs, density profile of the galactic 
halo). In addition, using the known cosmological abundance of dark matter in the Universe, 
we estimate the radiative WIMP production rate in the next generation of colliders (ILC) and 
compare their sensitivity to the WIMP mass with the XENON and GLAST projects. 

Our goal is not to perform an exhaustive analysis of the mass determination capacity of 
each detection mode. The GLAST discovery potential has been studied in much detail in ref.[9j 
whereas direct detection experiments are thoroughly treated in refs. [7] and [8]. Instead, we 
are mostly interested in examining at which point and under which circumstances the different 
detection techniques can achieve comparable precision, thus acting in a complementary way 
(an issue which has been commented upon in a great number of references but has never 
been treated as such). In the case of indirect detection, we will restrict ourself to W + W~ 
WIMP annihilation final states. It is well known, and is demonstrated in the Appendix, 
that different Standard Model particles contribute differently in the gamma-ray spectrum, 
the most prominent example being leptonic final states which tend to produce harder gamma- 
ray spectrae. Nevertheless, we should say in advance that at least in the ILC section we 
include an electron-positron final state, whereas the W + W~ spectrum can be considered 
quite representative of bosonic and hadronic ones. 

The paper is organised as follows. In section 2 we discuss the event rate and WIMP-nucleon 
scattering cross-section for a XENON-like experiment, in a microscopically model-independent 
approach. In Section 3 we carry out a similar analysis for the GLAST experiment, discussing in 
this case the WIMP annihilation cross section, and taking into account different halo profiles. 
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Section 4 is dedicated to the comparison between these two modes of detection. In Section 5 
we analyse the sensitivity that we can expect in such a model-independent framework for a 
linear collider. Finally, in Section 6 we carry out the comparison between the three detection 
modes. The conclusions are left for Section 7. 



2 Direct detection 
2.1 Differential event rate 

In spite of the experimental challenges, a number of efforts worldwide are actively pursuing 
to directly detect WIMPs with a variety of targets and approaches. Many direct dark matter 
detection experiments are now either operating or in preparation. All these experiments 
measure the number N of elastic collisions between WIMPs and target nuclei in a detector, 
per unit detector mass and per unit of time, as a function of the nuclear recoil energy E r . 
The detection rate in a detector depends on the density po ~ 0.3 GeV cm -3 and velocity 
distribution f{v x ) of WIMPs near the Earth. In general, the differential event rate per unit 
detector mass and per unit of time can be written as: 

dN _ a x - NP0 2 r° f(v x ) 

dE r ~ 2m 2 .™ 1 r) ' ~ * « ' 



F(E r y / J ~^dv x , (2.1) 



where the WIMP-nucleus cross section, o" x -tv, is related to the WIMP-nucleon cross section, 
<7 X _ P , by (J x -n = (J x _p(Am r / M T ) 2 , with M r = J^*^ the WIMP-nucleon reduced mass, 

n^r = m X +m N ^ ne WIMP-nucleus reduced mass, m x the WIMP mass, tun the nucleus mass, 
and A the atomic weight. F is the form factor. 

For the velocity distribution we take a simple Maxwellian halo 

f(^)d\ = ^^ 2 e-^ 2 d\, (2.2) 

where v x = 220 ± 20 km/s is the velocity of the Sun around the galactic center with its 
uncertainty, and we have neglected the motion of the Earth around the Sun. After integrating 
over the angular part in order to find the speed distribution we get: 

f(v x )dv x = J ^ 7 =e-M^ 2 dv x , (2.3) 

The integration over velocities is limited to those which can give place to a recoil energy E r , 
thus there is a minimal velocity given by v m i n (E r ) = "^f r ■ 

The effective interaction between the WIMP and a nucleus is given by the Woods-Saxon 
form factor 

m) = *MzM e -to.r , (2.4) 

qRi 

where the transferred momentum is q = \j2m^ E r , j% is a spherical, first-order Bessel func- 
tion, Ri = VR 2 -5 s 2 with R ~ 1.2 • A 1 / 3 fm, A is the mass number, and s — 1 fm. 

In FigHJwe show an example of a signal with a standard neutron background in a XENON- 
like (100 kg) experiment, after 3 years of data acquisition, as a function of the recoil energy. For 
a WIMP mass of 100 GeV and a WIMP-nucleon cross-section of 10~ 9 pb, such an experiment 
would reach a pretty large x 2 P er degree of freedom (Xred) > °f ^ ne or der of 60. 
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Figure 1: XENON event rate expectations for the case of Differential rate versus recoil energy of the nucleus for 

10~ 9 pb. The error bars shown are those expected for the 
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a WIMP mass m A 

XENON 100 kg experiment after 3 years of observation. The lower (blue) line is the background-only prediction. 
The x 2 P er degree of freedom (xj?cd) ls 59, giving a signal clearly distinguishable from the background. 



2.2 The statistical method 

In order to plot discrimination capacity regions, we use a method inspired by the treatment 
presented in [7]. For a WIMP with mass m r x eal and a WIMP-proton scattering cross-section 
v r x -p (in the case of indirect detection o~ x - v will, of course, be replaced by the total thermally 
averaged annihilation cross-section (crv)), we can calculate the theoretically expected number 
of events, say N t h, simply by integrating Eq. H2.in from E t h, the threshold energy which we 
consider as 4 keV, up to E sup , the maximal observable energy, which we take as equal to 30 
keV. Now, in a real-life experiment, one would expect the observed number of events to slightly 
deviate from this ideal value giving, say, Nexp events. In order to approach a more realistic 
situation, we perform our estimations not starting from Nth, but rather from Nexp, where the 
"experimental" number of events is drawn from a Poisson distribution with mean value Nth- 
Then, we can Monte-Carlo generate Nexp events according to Eq. (|2.ip and we obtain what 
could actually resemble to an experimental spectrum. Such a set of events, together with the 
corresponding value of Nexp will be in the following referred to as an "experiment". 

Then, for every point in the parameter space (m x ,a x - p ), we calculate the corresponding 
extended likelihood function which is given by: 



( N scan\N Exp Ne *p 

L = % i exp (~ N th an ) II a x- P ) (2-5) 



where 



f{h,m x ,a x -p) - E : (2.o) 



JZ up dN/dE(E;m x ,a 



x-pj 



is the normalized total event rate (signal+background) and N^ an is the theoretical number of 
events, expected from Eq. (|2.ip . for the given point of the parameter space. The normalization 
renders / a probability density function and, thus, suitable for use in a likelihood calculation. 
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The use of equation (12.51) presents the advantage that it takes into account the fact that 
the number of observed events in an experiment can, actually, deviate from the expected 
behaviour for several reasons. For the given experiment, say j, we scan over the (m x ,CT x _ p ) 
parameter space and calculate the value (m^ st '\ a^-p) that maximize the expression (12. 5p . 
This is the estimation for the j-th experiment. We then calculate the mean value of all the 
estimations and find which experiment's estimation was closest to this mean value. This 
experiment is considered to be the most representative of them all and is used to perform a 
final scan. Finally, from the likelihood distribution we obtain through this scan we can plot 
discrimination capacity regions. 

Direct Detection experiments present the advantage of quite well-controlled background. 
The additional ambiguity that arises in Indirect Detection and concerns uncertainties in the 
background will be dealt with in the relevant chapter. 

As a final remark on the statistical treatment we used, let us say that in order to be more 
precise, we would have to take into account (as is systematically done in [9]) the fact that the 
mass and cross-section precision are themselves random variables and should, consequently, be 
given with their relevant statistical variance. To do so, we would have to consider the actual 
distribution of estimators for all experiments. However, as discussed in the Introduction, 
such a treatment goes far beyond the scopes of this paper, where we are interested in a more 
qualitative comparison of different detection modes. In this respect, we keep the experiment 
which averages the properties of a larger set of experiments. Motivating this approach, our 
results are indeed in accordance with [9 J and |8j. 



2.3 The XENON experiment 

The XENON experiment at the Gran Sasso national laboratory aims at the direct detection of 
dark matter via its elastic scattering off xenon nuclei. It was deployed underground in March 
2006 and has been in continuous operation for a period of about one year. It allows the 
simultaneous measurement of direct scintillation in the liquid and of ionization, via propor- 
tional scintillation in the gas. In this way, XENON discriminates signal from background for 
a nuclear recoil energy as small as 4.5 keV. Currently a 10 kg detector is being used, but the 
final mass will be 1 ton of liquid xenon. In FigfJl we show the sensitivity curve for XenonlO 
(M = 10 kg) and Xenonl T (M = 1 ton) for T = 3 years of data acquisition, supposing zero 
backgrounds and a perfectly known velocity of the sun around the galactic center, fixed at 
220 km/sec. 

In our study, following Ref. [Tl] we will consider the energy range between 4 and 30 
keV and 3 years of data acquisition for a 100 kg XENON experiment. Such experimental 
conditions and time of exposure can be achieved after the 6 years of GLAST mission and 
justify the comparison between the two detection modes. 

In Fig|3]we show the ability of XENON to determine the mass and scattering cross-section 
for a 20 , 100 and 500 GeV WIMP. We can clearly see how sensitive the experiment is to light 
WIMPs: the precision can easily reach the percent level forQ m x < 50 GeV. Indeed, the recoil 
energy of the nucleus depends on the reduced mass (see Eq. (|2.ip ). For WIMPs much heavier 
than the nucleus mass (~ 100 GeV for Xenon), m r ~ uin, and is therefore independent of the 
WIMP mass. This is clearly reflected in the uncertainties at 68% and 95% CL in Fig|3]for a 
500 GeV WIMP. 

1 During the finalization of this work, Drees and Shan in Ref. [7] proposed that one can even increase such a 
precision with a combined analysis of two experiments of direct detection. 
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Figure 2: Spin-independent WIMP-nucleon cross-section versus WIMP mass for x 2 = 1,4 and M=10 kg and 1 
ton. 




Figure 3: Distribution of the maximum likelihood WIMP mass, m x , and cross-section, a x - p , for 3 years of 
exposure in a 100 kg XENON experiment, for m x = 20, 100, 500 GeV and <J X - P = pb. The inner (full) and 
outer (dashed) lines represent the 68% and 95% CL region respectively. The crosses denote the theoretical input 
parameters (a x - p , m x ). 
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2.4 Influence of astrophysical/background assumptions 

As mentioned before, a significant uncertainty (of the order of 8 — 10%) exists for the largely 
used values of the sun's circular velocity around the Galactic Center (GC), as well as the 
various background forms we could expect in direct detection experiments. In this respect, it 
would be interesting to examine how the previous results are altered in the case where vq is 
actually included in the fitting procedure, letting it vary within the given margin of error. 

As far as background events are concerned, it is quite difficult to perform a general study 
valid for every detector. Neutron backgrounds, which are in fact the most difficult to distin- 
guish from signal events, are usually taken to come from three sources (see also |13j): 

• Cosmic muon - induced neutrons, which are not in general considered to cause much 
nuisance. 

• Neutrons from the detector's surrounding rock. 

• Neutrons coming from contamination of the detector itself or surrounding 
materials. 

As we said, it is difficult to model in general neutron backgrounds, as they are mostly 
determined by the specific location in which every experiment is situated, as well as by the 
specific shielding configuration adopted by each collaboration. Two widely studied forms of 
neutron backgrounds are the case of a constant one, which seems to be quite well-motivated by 
an experimental point of view and can resemble to a heavy WIMP's signal, and an exponential 
one which apart from its theoretical motivation is also interesting as it gets to "mimic" (as 
pointed out in |8J) the actual signal spectrum for intermediate WIMP masses. In this respect, 
we studied the impact of these two forms of background: 

We consider firstly a constant background, with a value taken to be the same as the 
maximal WIMP signal in the first energy bin. Throughout this paper, when examining the 
impact of uncertainties on the mass determination accuracy, we will consider the case of a 
somehow "typical" in many theoretical frameworks case of a 100 GeV WIMP. 

Then, we introduce an exponential background of the form (^)^kg = Aex.p(— E/Ef,), 
where the slope of the exponential is fixed at E^ = 25 keV and the A factor is determined by 
demanding that the maximal values of the signal and the background be the same. The reason 
for this specific choice of parameters is that it is for these values that the signal spectrum has 
a significant resemblance to the background one, making it difficult to distinguish from one 
another. 

Our results are shown in Figs |4] and for the cases of inclusion of vq in the fitting procedure 
and non-zero backgrounds respectively. The deterioration of the expected accuracy is obvious, 
when we compare these plots to those of Fig|3l Especially for the case of large uncertainties 
in Vq (we let it vary in the whole margin-of-error region) and of inclusion of a background 
which is nearly degenerate with the signal, the expected precision is dramatically aggravated. 
This shows, among others, the extreme importance of a well-controlled environment and well- 
measured input parameters, other than the WIMP mass. 
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Figure 4: 68% and 95% CL regions for the XENON 100 kg experiment for a 100 GeV WIMP with a proton-WIMP 
scattering cross-section of 10 -8 pb in the case where uncertainties in the vq parameter are taken into account and, 
thus, included in the fitting procedure. 
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Figure 5: 68% and 95% CL regions for the XENON 100 kg experiment for a 100 GeV WIMP with a proton-WIMP 
scattering cross-section of 10 -8 pb including a constant neutron background (left) or an exponential one (right). The 
serious deterioration of accuracy in the second case is due to the fact that the background parameters where chosen 
in order to mimic the actual signal spectrum. 
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a (kpc) a (3 7 J(4 ■ 10~ 3 sr) 


NFW 
NFW C 
Moore et al. 
Moore c 


20 1 3 1 5.859 -10 2 
20 0.8 2.7 1.45 3.254 • 10 4 
28 1.5 3 1.5 2.574 • 10 4 
28 0.8 2.7 1.65 3.075 • 10 5 



Table 1: NFW and Moore et al. density profiles without and with adiabatic compression (NFW C and Moore c 
respectively) with the corresponding parameters, and values of J(Af2). 



3 Indirect detection 
3.1 Differential event rate 

The spectrum of gamma-rays generated in dark matter annihilations and coming from a 
direction forming an angle ip with respect to the galactic center is 

dN* i r 

^(E 7 ,i;)=J2 7 fBr i (av)—- I / p 2 dl , (3.7) 

unij oiun^ Jn ne f sight 

where the discrete sum is over all dark matter annihilation channels, dN^/dE^ is the dif- 
ferential gamma-ray yield, (av) is the annihilation cross-section averaged over its velocity 
distribution, Bri is the branching ratio of annihilation into final state "i" , and p is the dark 
matter density. The method followed in order to obtain the spectral function describing the 
standard model particle decay into 7-rays, is presented in the Appendix. 



It is customary to rewrite Eq. (|3.7|l introducing the dimensionless quantity J (which de- 
pends only on the dark matter distribution): 

J W = ^—( nor\l/ 3 ) / p\r{l^))dl. (3.8) 

8.5 kpc \0.3 GeV/cnrV J hne of sight 

After having averaged over a solid angle, Af2, the gamma-ray flux can now be expressed as 

$ 7 (£ 7 ) = 0.94 • 1(T 13 cnT 2 s- 1 GeV" 1 sr" 1 

. yffif r UH°gy J^Aii . (3.9) 

4" dE^ Vl0- 29 cm 3 s-V V m x J 

The value of J(AO)Af2 depends crucially on the dark matter distribution. The most common 
parametrization of the different profiles that have been proposed in the literature is 

P{ ) (r/Ro)f[l + (r/a)«]0»-7)/« ' 1 j 

where po is the local (solar neighborhood) halo density, a is a characteristic length, and Ro 
the distance from the Sun to the galactic center. As mentioned above, we will use po = 
0.3 GeV/cm 3 throughout the paper, but since this is just a scaling factor in the analysis, 
modifications to its value can be straightforwardly taken into account in the results. N-body 
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simulations suggest a cuspy inner region of dark matter halo with a distribution where 7 
generally lies in the range 1 (NFW profile [14]) to 1.5 (Moore et al. profile [E]), producing a 
profile with a behavior p(r) oc r -7 at small distances. Over a solid angle of 4 • 10~ 3 sr, such 
profiles can lead from J(Afi) ~ 5.859 • 10 2 to 2.574 • 10 4 . Moreover, if we take into account 
the baryon distribution in the Galaxy, we can predict even more cuspy profiles with 7 in the 
range 1.45 to 1.65 (J(AO) ~ 3.254 • 10 4 — 3.075 • 10 5 ) through the adiabatic compression 
process (see the study of Refs. [HUE]). We summarize the parameters used in our study and 
the values of J for each profile in Table [Q The values contained in the Table were calculated 
using the darkSUSY package for calculations of fluxes coming from SUSY DM candidates' 
annihilations. The calculation is obviously not altered whatever DM candidate is assumed 
and is therefore valid for an arbitrary WIMP as the ones examined here. 

It is worth noticing here that we are neglecting the effect of clumpyness, even though other 
studies showed that, depending upon assumptions on the clumps' distribution, in principle an 
enhancement of the flux by a factor of 2 to 10 is possible [18J. In this respect, the following 
predictions on the gamma-ray flux from the galactic center are conservative. 

3.2 Modeling the galactic center background. 

HESS [19J has measured the gamma-ray spectrum from the galactic center in the range of 
energy ~ [160 GeV-10 TeV]. The collaboration claims that the data are fitted by a power-law 

<P™ SS (E)=F E-: Y , (3.11) 

with a spectral index a = 2.21 ± 0.09 and F = (2.50 ± 0.21) • 10~ 8 m~ 2 s" 1 TeV" 1 . The 
data were taken during the second phase of measurements (July-August, 2003) with a \ 2 of 
0.6 per degree of freedom. Because of the constant slope power-law observed by HESS, 
it turns out possible but difficult to conciliate such a spectrum with a signal from dark 
matter annihilation [TTJEO] . Indeed, final particles (quarks, leptons or gauge bosons) produced 
through annihilations give rise to a spectrum with a continuously changing slope. Several 
astrophysical models have been proposed in order to match the HESS data [21]. In the 
present study we consider the astrophysical background for gamma-ray detection as the one 
extrapolated from the HESS data with a continuous power-law over the energy range of 
interest 1 - 300 GeV). As was recently underlined in Ref. [22J, the sensitivity of GLAST 
will be affected by the presence of such an astrophysical source. Note that the WIMP masses 
that we shall obtain in our parameter space < 1 TeV avoid any conflict with the observations 
of HESS. 

In addition, we have also taken into account the EGRET data [23] in our background 
at energies below 10 GeV (0^ RET (£')), as they can affect the sensitivity of the analysis. 
Indeed, the extrapolation of the gamma-ray fluxes measured by HESS down to energies as 
low as 1 GeV is likely to be an underestimation of the gamma-ray background in the galactic 
center, as EGRET measurements are one to two orders of magnitude higher than the HESS 
extrapolation. The EGRET point-source has been found [10J to be well fitted by a simple 
power-law with slope —2.2. We thus decided to take as background an interpolation between 
the HESS extrapolation and the EGRET data below 10 GeV to stay as conservative as possible 
in evaluating the gamma-ray background. 

Finally, we will consider the diffuse background of gamma rays in the region surrounding 
the galactic center. 

We will describe the spectrum of the background using the HESS observation from the 
Galactic Center Ridge [19], which can be described by 
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Figure 6: GLAST expectations for gamma-ray fluxes in the case of a WIMP mass m x = 180 GeV and cross 
section (av) = 3 • 1CP 26 cm s . A NFW halo profile has been adopted. The error bars shown are those projected 
for the GLAST experiment after a six-year mission run, assuming the galactic center will be within its field-of-view 
50% of the time. The lower (blue) line is the background-only prediction. The xled ls 132, giving a signal clearly 
distinguishable from the background. 



^bkgC^O = 11 ■ 10- 4 ^GeV 9GeV_lcm ~ 2s " lsr_1 ' ( 3 - 12 ) 

In our analysis, we will consider the inner 2° x 2° field of view (Afi = 4- 10~ 3 sr) and the energy 
region between 1 and 300 GeV. During the completion of our work, the authors of Ref. [10] 
gave a more detailed and sophisticated statistical analysis of the diffuse background, adopting 
an overall normalization around the galactic center and taking into account a statistical spread 
function as GLAST would be able to probe in this region. However, we have checked that our 
results are not significantly modified and we recover similar results concerning the prospects 
of GLAST. 

3.3 The GLAST experiment 

The space-based gamma-ray telescope GLAST [12] was launched in May 2008 for a five-year 
mission. It will perform an all-sky survey covering a large energy range (« 1 - 300 GeV). 
With an effective area and angular resolution of the order of 10 cm 2 and 0.1° (AJ7 ~ 10~ 5 
sr) respectively, GLAST will be able to point and analyze the inner center of the Milky Way 
(~ 7 pc). Concerning the statistical method, we have used an analysis similar to the one 
considered in the case of direct detection in section l2~2l with a six- year mission run, assuming 
the galactic center will be within its field-of-view 50% of the time [ElESJ. In FigEJwe show 
the ability of GLAST to identify a signal from dark matter annihilation for a WIMP mass of 
180 GeV. The error bars shown are projected assuming Gaussian statistic, and we adopt the 
background described above including Poisson noise. In the following, we will concentrate on 
a process which gives 100% annihilation to WW. We have checked that the dependence on 
the final state does not influence significantly the general results of the study, except in the 
case of leptonic final states. This will be studied in a specific case in Section 5 dedicated to 
the ILC experiment. 

In FiglT] we show the ability of GLAST to determine the mass and annihilation cross- 
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Figure 7: Distribution of the maximum likelihood WIMP mass, m x , and annihilation cross-section, (av), after 6 
years of observation (50% of time exposure) of the galactic center with GLAST, with the hypothesis of a NFW halo 
profile, for m x = 50, 100, 500 GeV and (av) = 3 • 10~ cmV 1 . The inner (full) and outer (dashed) lines represent 
the 68% and 95% CL region respectively. The crosses denote the theoretical input parameters ((av), m x ). 



section for a 50, 100 and 500 GeV WIMP. Again, we can see that the experiment is more 
sensitive to light WIMPs: the precision can easily reach the percent level for GLAST for 
m x < 50 GeV. The gamma-ray spectrum will give more precise measurements if the mass 
of the WIMP lies within the GLAST sensitivity range. Indeed, the shape of the spectrum 
will be easily reconstructed above the HESS/EGRET and diffuse background if the endpoint 
of the annihilation spectrum lies within the energy range reachable by GLAST [0.1 — 300 
GeV]. Furthermore, we have studied the influence of the variation of the inner slope of the 
halo profile on the resolution of the WIMP mass. In addition to the NFW profile, we have 
considered some NFW-like profiles, allowing the 7 parameter in Eq. ()3.10|) to vary from its 
original value by 10%. This is shown in Fig|8j where in addition to the NFW halo profile 
(7 = 1) we also study profiles with 7 = 0.9, 1.1. As expected, the larger the 7 is, the more 
enhanced the galactic gamma ray flux becomes, and the better the WIMP mass resolution 
turns out to be. It is worth noticing here that in the case of a compressed NFW profile 
(7 ~ 1-45), the precision of GLAST increases by two orders of magnitude. 

3.4 Varying background parameters 

Throughout the previous (and the following) analysis, we have considered a perfectly known 
background for both detection modes. Whereas this is a rather reasonable approximation for 
the case of direct detection, it is less obvious for the indirect one. As it has been pointed 
out (see, for example, [SUES]), the uncertainties entering the calculation of the backgrounds 
coming from the galactic center region can considerably affect the results of any analysis. More 
concretely, and especially for small WIMP masses and low energies (where the performance 
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Figure 8: NFW-like halo profile with 7 = 0.9, 1 (NFW) and 1.1 at 95% confidence level. 



of both direct and indirect detection is maximal), the main contribution in background comes 
from the EGRET source aforementioned. However, both the overall normalisation and the 
spectral index characterising this source's spectrum contain uncertainties, and it has been 
widely discussed in the literature that a dark matter annihilation signal could account for 
some part of the EGRET source signal. An interesting point would be to include the overall 
background normalization as well as the spectral index in the fitting procedure. In Fig|9] 
we show, for the sake of comparison, the result of a fitting procedure, where we also fit 
the background normalization while simultaneously considering signals and backgrounds with 
poissonian fluctuations. The original spectrum is taken to be the full EGRET source plus the 
flux produced by a 100 GeV WIMP annihilating in a NFW halo. One could imagine discarding 
low-energy data which contain a maximal background contamination. This, however, would 
significantly reduce the statistics and the corresponding precision, since a major part of the 
signal would be discarded. The fact that the inclusion of an uncertainty in the background 
normalization (i.e. its inclusion in the statistical treatment) does not have a major impact on 
the results can be explained from the fact that throughout this work we have used eq. (|2.5p 
in our statistical analysis, which already introduces a deviation from the ideally expected 
situation. In this respect, our results are quite conservative. 

In the same way, in FigJT0]we show the corresponding results where this time the spectral 
index is included in the fitting procedure instead. The spectral index is left to vary in the 
region [2.1,2.4], which we find to be a quite reasonable one as we verified that all over this 
region we obtain reasonable fits of the EGRET data. 

It is interesting to note that the variation of the background's spectral index seems to have 
a larger impact on the precision that could be achieved, with respect to the corresponding 
case of the background's overall normalisation. 

This is somehow logical, first of all since by definition the background depends linearly 
on the normalisation factor, but exponentially on the spectral index of the EGRET point 
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Figure 9: 68% and 95% CL regions for a statistical treatment with the overall background normalization included 
in the fitting procedure, m x = lOOGeV and (av) = 3 ■ 10 -26 cm 3 s _1 . 
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Figure 10: 68% and 95% CL regions for the case where the uncertainties in the EGRET point source spectral 
index are included in the fitting procedure, m x = lOOGeV and (av) — 3 • 10~ 26 cm 3 s _1 . An NFW halo profile has 
been assumed. 
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source. So, modifications in the latter bring along a much more drastic modification of the 
background signal itsself. Furthermore, variations of the overall normalisation have just the 
influence of "burying" the signal a little more or a little less in a background which is already 
quite elevated. On the contrary, by varying the spectral index we actually change the shape 
of the spectrum. This brings along a more important uncertainty, since we could imagine 
much more numerous configurations in the (spectral index, cross-section, mass) space that 
could satisfy selection criteria. 

3.5 Impact of Final States 

Again, throughout this paper, we have considered a pure WIMP annihilation into a W + W~ 
final state. This is an assumption which is made to simplify the overall treatment, but which 
at the same time somehow restricts the generality of our results. In this paragraph, we are 
interested in examining what could be the impact of variations in the final state on the WIMP 
mass determination capacity. Annihlation into ZZ pairs is not expected to seriously modify 
the results, what would be more interesting would be to see what happens when we consider 
(light or heavy) quark pairs and/or leptons as WIMP annihilation products. 

The spectrae of SM particles' annihilations into gamma-rays can be found in the Appendix. 
As explained there, the only leptonic final state we consider is the r + r _ one, since annihila- 
tion into has a relatively small contribution to the annihilation gamma-ray spectrum, 
whereas e + e~ pairs contribute through other processes, the examination of which exceeds the 
purposes of the present treatment. We should, nevertheless, note here that we do not take 
into account the effects of leptonic final state radiation which can indeed become important, 
especially in the case of Kaluza-Klein dark matter and in energy ranges lying near the WIMP 
mass. The effect of such processes has been discussed in detail in refs.[27j for the case of 
a generic WIMP and [28J for the special case of KK dark matter. Obviously, this omission 
somehow retricts the generality of our results as far as the impact of final states are concerned. 

To this goal, we performed two kinds of tests: The first one consists only of modifying the 
annihilation products, considering a perfectly known final state (meaning that the Branching 
Ratios are not included in the statistical treatment). Our results can be seen in FigJTT]for the 
case of pure bb, qq and t + t~ final states and a 100 GeV WIMP. 

It is interesting to notice the relative amelioration of the mass resolution with respect to 
the pure W + W~ final state. A possible explanation could be that, as can clearly be seen in 
gamma-ray yields presented in the Appendix, fermionic final states tend to give more hard 
spectrae with respect to bosonic ones (the extreme case being leptonic ones), rendering the 
spectrum more easily distinguishable from the background. The hardest spectrum is given by 
the t + t~ final state. Nevertheless, in this case the characteristic spectral form is somewhat 
compensated from the reduced statistics of the signal. This is not the case for annihilation 
into quarks, where the characteristic spectral form, although less obvious than in the leptonic 
case, is nevertheless combined with an important enhancement in the signal. 

As a second test, we consider a mixed final state and fit the BRs themselves along with 
the annihilation cross-section and the WIMP mass. Our results can be seen in FiglT2l where 
we have taken a 100 GeV WIMP annihilating into a final state consisting of 70% l^+W - and 
30% t + t~ . The sum of the two branching fractions is obviously equal to 1, so we only need 
to include one further parameter in the statistical analysis. 

In this case, we can clearly see that the mass resolution deteriorates with respect to the 
case where a perfectly known final state is considered. A possible explanation could be that a 
mixed lepton - gauge boson final state yields a spectrum which presents neither the augmented 
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Figure 11: 95% CL regions for a 100 GeV WIMP and different final states. The total thermally averaged 
annihilation cross-section has been taken to be (<rv) — 3 • 10~ 26 cm 3 s _1 . A NFW halo profile has been assumed. 

statistics of pure annihilation into gauge bosons (the gamma-ray yield of leptons is significantly 
inferior to the one of gauge bosons) nor the characteristic hard spectral form of annihilation 
into t + t~ pairs. 

4 Direct versus indirect detection experiments 

It is interesting to remark from Figsj3]and[7l how two completely different means of observa- 
tion, with completely different signal/background physics, are in fact competitive (and hence 
complementary) in the search for the dark matter. 

In Fig(l3]we compare the precision level for both experiments as a function of the WIMP 
mass, for different values of the spin-independent cross-section (10~ 7 , 10~ 8 and 10 -9 pb) 
and for different halo profiles. For this treatment we minimize the impact of uncertainties 
discussed in the previous paragraph, as we are mostly interested in examining the a priori, in 
some sense "intrinsic" sensitivity of the two detection techniques. For example, in the case of 
direct detection, the necessity for minimisation of background noises and control of the noise 
sources has been repeatedly stressed out. As for uncertainties entering the velocity distribution 
of WIMPs in the solar neighbourhood (or, why not, the form factor's functional form), these 
can, in principle be minimized by measurements exterior to the experiments themselves. 
The same holds for uncertainties in the case of indirect detection. As a small example, 
Fermi-GLAST's overall sky survey capacity is hoped to contribute in the minimisation of 
uncertainties in non-DM annihilation sources, whereas other observations in different energy 
regions can also contribute in this direction. In this respect, for our comparative results, we 
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Figure 12: 68% and 95% CL regions for a 100 GeV WIMP with (err;) = 3- 10~ 26 cm 3 s" 1 considering a mixed final 
state consisting of 70% W + W~ and 30% t + t - and including the Branching Ratios in the statistical treatment. 



remove the extra parameters from the statistical treatment. We see that at 95% of confidence 
level GLAST, after 3 years of exposure (6 years of taking data at 50% of time exposure), will 
have an equivalent sensitivity to a 100 kg XENON-like experiment after 3 years of running 
if o" x _p < 10~ 8 pb, independently of the WIMP mass. The indirect detection by GLAST will 
always be able to give an upper bound on the WIMP mass for m x ~ 100 GeV, whereas a 
XENON-like 100 kg experiment would only give a lower bound value if o x - v < 10~ 9 pb. In 
all cases, the lower bounds given by GLAST for a NFW halo profile are similar to the ones 
given by a XENON-like 100 kg experiment for any WIMP mass if <J X - P < 10 -8 pb. 

To compare the uncertainties on the WIMP mass expected from direct and indirect detec- 
tion modes, we plotted in FigfT4l as a function of the WIMP mass for o x - v = 10 _8 pb and 
a NFW halo profile. One can clearly see in the figure that GLAST will be competitive with 
XENON 100 kg to measure the WIMP mass in the case of a NFW halo profile for a x - p < 10~ 8 
pb. 

5 WIMPs at present and future colliders 

Among the most important sources of information concerning WIMP dark matter are, obvi- 
ously, collider experiments, both present, such as the Tevatron or the Large Hadron Collider 
(LHC) and future, such as the International Linear Collider (ILC). In fact, there is a quite 
general agreement on the fact that despite the significant progress in astroparticle physics 
experiments, which manage to impose more and more constraints on various models, collider 
experiments remain an irreplaceable source of information for particle physics. It is quite 
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Figure 13: Comparison between 100 kg XENON-like experiment and GLAST sensitivity in the case of different 
halo profiles, at 95% of confidence level, for several WIMP masses (50, 100 and 500 GeV) and WIMP-nucleon 
cross-sections (10 -7 , 10 -8 and 10 -9 pb). In each panel the cross denotes the input parameters. 
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natural thus, to examine the potential of colliders to constrain WIMP properties. We will ex- 
amine the extent at which astroparticle and collider experiments become competitive, trying 
at the same time to stay as model-independent as possible. 

This last point is, in fact, the major difficulty in treating collider experiments to extract 
astroparticle data: most studies performed for new physics at colliders are very strongly model 
dependent. This is almost unavoidable for the case of the LHC, due to the hadronic nature 
of the colliding particles. The large uncertainties on the parton distribution functions (and, 
hence, on the initial energy of the colliding particles/partons) render it extremely difficult (in 
fact, practically impossible) to look beyond the transverse plane. This fact obviously limits 
-up to a certain point- the precision that could be obtained with respect to, for example, an 
e + e~ collider. As a result, it is quite difficult to make predictions in a model-independent 
way, since a whole set of parameters must be taken into account in order to perform concrete 
predictions. The cruciallity of these uncertainties will become clearer in the following. 



5.1 The Approach 

Recently, an approach was proposed in references [29j [30] which allows to actually perform a 
model-independent study of WIMP properties at lepton colliders (such as the ILC project). 
The goal we pursue is to extract constraints which are as stringent as possible for a generic 
dark matter candidate. A generic WIMP can annihilate into pairs of standard model particles: 

X + x-^Xi + Xi . (5.13) 

However, the procedure taking place in a collider is the opposite one, having only one species 
of particles in the initial state. The idea proposed in Ref. [29] is to correlate the two processes 
in some way. This can be done through the so-called "detailed balancing" equation, which 
rGctds * 

a(x + X^X l + Xj) _ v 2 x (2S x + l) 2 

a(X l + X t ^ X + x) v 2 x (2S x + ir ' { ' ' 

where Vi and Si are respectively the velocity and the spin of the particle i. The cross-section 
a (xX ~^ XiXi) is only averaged over spins. 
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The total thermally averaged WIMP annihilation cross-section can be expanded as 

oo 

E^ (JVJ ' ( 5 - 15 ) 



OiV = 

J=0 



where J is the angular momentum of each annihilation wave. Now, for low velocities, the 
lowest-order non-vanishing term in the last equation will be dominant. So, we can express 
the total annihilation cross-section as a sum of the partial ones over all possible final states 
for the dominant partial wave Jq in each final state: 



-£^ /0) - (5-16) 



Next, we can define the "annihilation fraction" Ki into the standard model particle pair 

X{ — X{\ 

a (Jo) 

Ki = -!— . (5.17) 
Can 

By combining Eqs. (|5.14|) and (15.17P we can obtain the following expression for the WIMP 
pair-production cross-section: 



oiX X - 2vl a ( 2g x + I) 2 L 1/2+J ° 

a(X lXl - 2 X ) = 2 ] ^ {2Sx + l? II " — I 



(5.18) 



Now, a few remarks should be made about the validity of this formula: 

• Equation (15, 18ft is valid for WIMP pair-production taking place at center-of-mass ener- 
gies just above the pair-production threshold. 

• The detailed balancing equation is valid if and only if the process under consideration 
is characterized by time-reversal and parity invariance. It is well known that weak 
interactions violate both of them, up to some degree, which we ignore in this treatment. 

A process of the form XiXi — ► \X is n °t visible in a collider, since WIMPs only manifest 
themselves as missing energy. At least one detectable particle is required for the event to 
pass the triggers. An additional photon from initial state radiation (ISR) is required to be 
recorded on tape: XiX,- t — > XXI- We can correlate the WIMP pair-production process to the 
radiative WIMP pair-production for photons which are either soft or collinear with respect to 
the colliding beams. In this case, the two processes are related through |29j : 

da{e+ Tr 2 a + 7) « Fix, cos 0)<7(e+e- - 2 X ) , (5.19) 
axa cos v 

where x = 2E^/y/s, 9 is the angle between the photon direction and the direction of the 
incoming electron beam, a is the WIMP pair-production cross-section produced at the reduced 
center of mass energy s = (1 — x)s, and T is defined as: 

^ al + (l-x) 2 1 , 

F(x,cos9) = i 5.20 

7r x sin 6 

Now, by combining Eqs. (15, 19p and (|5.18h we get the master equation: 

jg-i&o- - 2 x+ 7) - + 1 2 1 1 _ 4M|\ UM " 

dxacosU 167r x sin \ (1 — x)s I 

(5.21) 
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The problem is that very collinear photons fall outside the reach of any detector, due to 
practical limitations in the coverage of the volume around the beam pipe. Also, typically, 
lower cuts are included in the detected transverse momentum of photons, pt = E 7 sm9, in 
order to avoid excessive background signals at low energies. So, if we are to use this approach, 
we have to examine its validity outside the soft /collinear region. The accuracy of the collinear 
approximation for hard photons at all angles has been discussed in the original paper [29], 
with the conclusion that the approach works quite well. 

However, an important point should be taken into account here. Prom the previous dis- 
cussion on the validity of the method, we have to impose specific kinematical cuts on the 
detected photons. We consider the following conditions: 

• We demand an overall condition sin6 > 0.1 and px > 7.5 GeV in order to assure the 
detectability of the photons. 

• In order to assure the fact that any photon under examination corresponds to non- 
relativistic WIMPs, we demand v\ < 1/2. This gives a lower kinematical cut, along 
with an upper cut corresponding just to the endpoint of the photon spectrum: 

T^—^j^T^—^)- (5 22) 

These conditions present a flaw: the energy limits depend on the mass we wish to constrain. 
On the other hand, for the reasons explained before, we cannot treat the signals without 
imposing such kinds of cuts, if we do not want either to abuse the method or stick to heavy 
WIMPs (which, for kinematical reasons, cannot be relativistic) . The only way to evade this 
problem is to suppose that other dark matter detection experiments (or, eventually, the LHC 
in the framework of specific models) will have already provided us with some sort of limits on 
the WIMP mass. In this case, having an idea of the region in which the WIMP mass falls, 
we can also estimate the cuts that will safely keep us outside the relativistic region and only 
consider photons within this region. 

The main source of background events is the standard model radiative neutrino production, 
e + e~ — > vvj. Apart from these background events, various models predict additional signals 
of the form "7 + missing energy", one of the most well-known examples being radiative 
sneutrino production [311 132| . predicted in the framework of several supersymmetric models. 
In the spirit of staying as model-independent as possible, we will ignore all possible beyond 
standard model processes. 



5.2 Basic Results 
5.2.1 Non-polarized beams 

We place ourselves in the framework of the ILC project with a center-of-mass energy of 
y/s = 500 GeV and an integrated luminosity of 500 fb^ 1 . In order to estimate the background 
events, we used the CalcHEP code (33j [34] to generate 1.242.500 e + e~ — > i/vj events, 
corresponding to the aforementioned conditions. The total radiative neutrino production 
background can be seen in FigHH The peak at E 1 = yfs/2 • (1 - Mf/s) ~ 241.7 GeV 
corresponds to the radiative returns to the Z resonance. 

We generate a predicted "observable" spectrum for given values of the WIMP mass and the 
annihilation fraction. During this study, we do not proceed to a (more realistic) full detector 
simulation, as done for example in Ref. [35], but stick to prediction levels in order to perform 
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Figure 15: Radiative neutrino production background e + e — > vw) for the ILC, for an unpolarized initial state. 



as thorough a scan as possible in the (m x ,K e ) parameter space and to have a picture of the 
"a priori" potential of the method. 

Figure [16] shows the predicted ability of the ILC to determine WIMP masses and annihi- 
lation fractions for four points in the (m x , n e ) parameter space for a 68% and 95% Confidence 
Level. These results concern WIMPs with spin S x = 1/2 and an angular momentum Jo = 1 
which corresponds to an annihilation cross-section a an = 7 pb [29]. As can be seen, we are 
able to constrain quite significantly the WIMP mass (20% — 40% precision), while constraints 
on K e are weaker. 

Figure [T71 shows the relative error (Am x /m x ) for the mass reconstruction as a function 
of m x , for n e = 0.3 and a 95% confidence level. The solid line corresponds to the proper 
treatment including kinematical cuts. For indicative reasons, we also show the abused results 
obtained if we do not impose kinematical cuts on the photon energy (dashed line). The 
amelioration of the method's efficiency is obvious, although this is after all a false fact, since 
we include regions in which the approach is not valid. Above m x ~ 175 GeV the two lines 
become identical, since the WIMPs cannot be relativistic. The capacity of the method peaks 
significantly for masses of the order of 175 GeV because around this value we reach an optimal 
combination of phase space volume and available spectrum that passes the kinematical cuts 
and can, hence, be used for the calculation of the relevant x 2 , whereas as we move away from 
this value the accuracy tends to fall. 

Let us make a final remark on the possibility of adopting a similar approach in the case 
of the LHC. As we argued before, the large uncertainties entering the parton distribution 
functions and, hence, the large uncertainty in the collision energy, affect significantly the 
precision of the whole procedure (which is, already, based on approximations). Formally, in 
order to perform such a study for the LHC, the computed cross-sections must be convoluted 
with the proton form factors. As an additional element, the photon background in the LHC is 
expected to be much greater than in the ILC. The possibility of determining WIMP properties 
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Figure 16: Distribution of the maximum likelihood WIMP mass and annihilation fraction for the ILC in the 
(rn x ,K e ) plane, for n e — 0.3 and m x = 125,150,175 and 200 GeV. The inner (full lines) and outer (dashed lines) 
represent the 68% and 95% CL region respectively. 
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Figure 17: Relative error in a generic WIMP mass determination, for n e = 0.3 and at a 95% confidence level. 
The solid line corresponds to the results obtained after imposing the proper kinematical cuts, whereas the dashed 
line to the case where we do not take these limits in consideration. 
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Figure 18: Relative error in a generic WIMP mass determination, for three cases of beam polarization, including 
all proper kinematical cuts. 



through a model-independent method at the LHC has been addressed to in Ref. [36], where 
the authors conclude that WIMP detection will be extremely difficult, if even possible. 

5.2.2 Polarized beams 

The reach of the ILC can be further increased by polarizing the beams. For polarized beams, 
the signal cannot be fully characterized by K e ; instead, four independent annihilation fractions 
are needed, corresponding to the four possible e + e~ helicity configurations. 
To apply Eq. (|5.21|) to this case, we make the replacement: 

Ke - \(l + P-) [(l + P + )n{e^e L + ) + {l-P + )K(e R el)] 

+ \{1 - P_) [(1 + P+) K (e L _e L + ) + (1 - P+) «(e£e?)] , (5.23) 

where P± are the polarizations of the positron and the electron beams. As in ref [29j[35], let us 
assume that the WIMP couplings to electrons conserve both helicity and parity: K(e^e+) = 
K^tef) = 2n e and K(e^e^) = K{e L _e\) = 0. 

In Fig{l8] we show the relative error for the mass reconstruction for K e = 0.3 and 95% 
confidence level, for the unpolarized scenario and for two different polarizations: (P_,P + ) = 
(0.8,0) and (0.8,0.6). 

6 Complementarity 

In Fig(19]we compare the precision levels for direct and indirect detection experiments, along 
with the corresponding results of the method we followed for the ILC for two cases of WIMPs 
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m x 


XENON 


GLAST 


ILC 


50 GeV 
100 GeV 
175 GeV 
500 GeV 


-5/ + 7 GeV 
-19/ + 75 GeV 
-65/ GeV 


±12 GeV 
-50/ + 60 GeV 
-125 GeV 


-40/ + 20 GeV 
-20/ + 15 GeV 



Table 2: Precision on a WIMP mass expected from the different experiments at a 95% CL after 3 years of exposure, 
<7 X _ P = 1CP 8 pb a NFW profile and a 500 GeV unpolarized linear collider with an integrated luminosity of 500fb~ 1 

masses, m x = 100 GeV and 175 GeV, and n e = 0.3. We plot the results in the {m x , K e ) plane. 
This is done as the K e parameter entering the ILC treatment presented before is, in fact, the 
same parameter as the corresponding branching ratio Bri = appearing in eq. (13. 7p for 
i = e. 



The blue-dotted line corresponds to a 100 kg XENON- like experiment, where the WIMP- 
nucleus cross-section has been assumed to be 10~ 8 pb. The green-dashed line depicts the 
results for a GLAST-like experiment assuming a NFW halo profile. The total annihilation 
cross-section into standard model particles has been taken to be (av) = 3- 10 -26 cm 3 s _1 . The 
red-plain line represents our results for an ILC-like collider, with non-polarized beams. All 
the results are plotted for a 95% confidence level. 

We can see that for different regions of the WIMP mass, the three kinds of experiments 
that we have used as prototypes can act in a highly complementary way. For example, for the 
case of a 100 GeV WIMP, indirect detection or an ILC-like experiment alone can provide us 
with limited precision both for the WIMP mass (of the order of 60%) and the K e parameter 
(where the results are even worse). Combined measurements can dramatically increase the 
precision, reaching an accuracy of 25% in mass. If we additionally include direct detection 
measurement, we can further increase the precision. 

For the case of a 175 GeV WIMP, a point where the unpolarized ILC sensitivity peaks, 
we see that the dominant information comes from this source. Nevertheless, even if we only 
combine direct and indirect detection experiments, we see that we can, in fact, acquire non- 
negligible constraints on the dark matter candidate mass. 

To summarize the analysis, we show in Table[2]the precision expected for several interesting 
dark matter masses. Whereas a light WIMP (50 GeV) can be reached by both types of dark 
matter experiments with a relatively high level of precision, our analysis fails in the ILC case 
because of the relativistic nature of the WIMP. On the contrary, the ILC would be particularly 
efficient to discover and measure a WIMP with a mass of about 175 GeV. Concerning a 500 
GeV WIMP, which is kinematically unreachable at the linear collider, it would be difficult to 
be observed by GLAST or XENON. Only a lower bound could be determined experimentally. 

7 Conclusions 

A Weakly Interacting Massive Particle (WIMP), with mass lying from the GeV to the TeV 
scale, is one of the preferred candidates for the dark matter of the Universe. 

We have discussed the possibility of identifying WIMP properties in a model-independent 
way. For that we have considered direct and indirect searches, and in particular the interest- 
ing cases of a XENON-like 100 kg experiment and the GLAST satellite. We have shown that 
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Figure 19: Comparison between a 100 kg XENON-like experiments (dotted line) with er x _ p = 10~ 8 pb, GLAST 
(dashed line) in the case of an NFW halo profile with (av) = 3 • 10~ 26 cm 3 s _1 , and unpolarized ILC sensitivity 
(solid line) at 95% of confidence level, for different WIMP masses m x — 100 and 175 GeV, and n e = 0.3. 
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whereas direct detection experiments will probe efficiently light WIMPs given a positive de- 
tection (at the 10% level for m x < 50 GeV), GLAST will be able to confirm and even increase 
the precision in the case of NF W profile, for a WIMP-nucleon cross-section (J x - V < 10~ 8 pb. 

Moreover, both XENON and GLAST are complementary with a future ILC project, and 
the measurements from the three experiments will be able to increase significantly the precision 
that we can reach on the mass of the WIMP. 
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Appendix 



In this Appendix we present the method followed in order to obtain the functions describing the 
standard model particle decay into 7-rays. In order to determine these spectral functions, we 
generated 300000 events of standard model particles decaying (directly or through secondary 
decays) into 7-rays using the PYTHIA [37J package, taking care in order to include all possible 
decay channels. Following the method of Ref. [38], we fitted the resulting spectra through 
functions of the form: 

dN l 

= exp[F i (ln(x))], (7.24) 

where i represents the i-th WIMP annihilation channel, i = WW, ZZ, etc; x = E^/m x with 
m x being the WIMP mass and F are seventh-order polynomial functions which were found 
to be the following: 

WW(x) = -7.72088528 - 8.30185509 x - 3.28835893 x 2 - 1.12793422 x 3 

- 0.266923457 x 4 - 0.0393805951 x 5 - 0.00324965152 x 6 - 0.000113626003 x\ 
ZZ{x) = -7.67132139 - 7.22257853 x - 2.0053556 x 2 - 0.446706623 x 3 

- 0.0674006343 x 4 - 0.00639245566 x 5 - 0.000372241746 x 6 - 1.08050617 • 10~ 5 x 7 , 
bb(x) = -11.4735403 - 17.4537277 x - 11.5219269 x 2 - 5.1085887 x 3 

- 1.36697042x 4 - 0.211365134x 5 - 0.0174275134 x 6 - 0.000594830839 x 7 , 
uu(x) = -4.56073856 - 8.13061428 x - 4.98080492 x 2 - 2.23044157 x 3 

- 0.619205713 x 4 - 0.100954451 x 5 - 0.00879980996 x 6 - 0.00031573695 x 7 , 
dd{x) = -4.77311611 - 10.6317139 x - 8.33119583 x 2 - 4.35085535 x 3 

- 1.33376908 x 4 - 0.232659817 x 5 - 0.0213230457 x 6 - 0.000796017819 x 7 , 
t -t+(x) = -5.64725113 - 10.8949451 x - 7.84473181 x 2 - 3.50611639 x 3 

- 0.942047119x 4 - 0.14691925x 5 - 0.0122521566 x 6 - 0.000422848301 x 7 . 

The case of WIMP annihilation into pairs has a relatively small decay contribution, 

to the photon spectrum, coming from the fj, — > e~v e y il pj channel, which has a small branching 
ratio. e + e~ pair production contributes to the gamma-ray spectrum through different (not 
decay) processes, mainly inverse Compton scattering and synchrotron radiation. These contri- 
butions depend crucially on the assumptions made concerning the intergalactic medium and 
will not be analysed here (For a relevant treatment see, e.g., section V of [39J ) . This means, 
practically, that the e + e~ and fi + pT spectral functions are set equal to zero. A graphical 
representation of these functions can be seen in Figl20l 

These functions can afterward be used in order to generate any gamma-ray flux according 
to eq. ([377) 

As we can see, all contributions are quite similar, apart from the r~r + channel which 
has a characteristic hard form. Nevertheless, at high energies, the form of all contributions 
becomes almost identical. 
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Figure 20: Separate contributions from standard model particles decaying into 7-rays for m x — 100 and 500 GeV. 
The PYTHIA result points have been suppressed for the sake of clarity. 
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